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1. Introduction 

In the quest to explore physics beyond the Standard Model (SM), certain new 
proposals have held considerable interest of the community. Amongst these is the 
proposal by Randall and SundrumP] wherein the SM fields live on one of the two 3-branes 
which are end of the world branes of a five dimensional spacetime. This scenario can be 
embedded in higher dimensional spacetime as well. One of the attractive features of this 
proposal is the resolution of the hierarchy problem between the Planck scale and the elec- 
troweak scale. This is achieved by choosing the geometry of the embedding spacetime with 
an exponential warp factor for the four dimensional spacetime which is the worldvolume of 
the 3-brane. This exponential warp factor produces a difference in the mass scales between 
the two end of the world 3-branes. 

There are several variants of the Randall- Sundrum (RS) model involving multiple 
branes intersecting brane configurations and supersymmetry[0-p. Many of these 
models including the RS model have atleast one negative tension brane. This is demanded 
by the requirement of charge neutrality in the compact space. For example, in the original 
RS model, the visible world lives on the negative tension brane {i.e., the one with the 
induced cosmological constant Ams < 0). 

Any proposal for physics beyond the SM has to pass very stringent tests laid down 
by the SM itself as well as by the cosmological observations that have been made till date. 
Brane world cosmology has been studied in refs.[|T0|-|T7|. These papers conclude that the 
rate of expansion of the universe in the RS model is different from that of the familiar 
Friedman-Robertson- Walker cosmology in four dimensions. In particular, in RS models, 
the Hubble parameter is proportional to the energy density p, i.e., H ^ p. In contrast, in 
our universe, the Hubble parameter seems to behave as if ~ ^/p. It is, however, possible 
to circumvent this problem if one realizes that the total energy density p is the sum of the 
vacuum energy density, i.e., the brane tension and the matter energy density pm which 
lives on the brane. In this case, dependence of the Hubble parameter on the energy density 
is if ~ K^isPm- This resolution brings up a new problem though. If the RS model is to 
solve the hierarchy problem, the visible brane needs to have a negative tension. It turns out 
that the Hubble parameter on such a brane could be real only if the matter energy density 
is negative. However, this would lead to antigravity in the brane world, a consequence 
clearly at variance with experimental observations. 

Another issue which was not dealt with in the original RS proposal was that of 
the stability of end of the world brane model. This was resolved by Goldberger and 
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Wise (GW)[jT^^. They introduced a bulk scalar field into the model and showed that its 
(minimal) coupling to the bulk gravity stabilizes the two end of the world 3-branes at a 
critical distance r^. For all parameters of o(l) this can generate a warp factor of the order 
10^^ GeV thereby relating the Planck Scale to the SM scale. 

Recently, Kogan et al.[^^ proposed a modified version of the RS model which contains 
three 3-branes with both the end of the world branes (one at y = and the other at y = L2) 
being of positive tension and the third (moving) brane with negative tension (at y = Li). 
We will call this the three 3-brane model. The visible world lives on one of the end of the 
world branes (at y = L2). This construction automatically ensures that the cosmological 
constant on our brane is positive and thereby avoids the problem of the Hubble parameter. 
This model has an exponential warp factor between either of the positive tension branes 
and the negative tension brane. The phenomenological consequences are rather striking 
and have been studied in refs.p,|2^. 



Stability of the Kogan et al. [|2^ proposal is quite crucial. This is because the moving 
brane has negative tension and we have exponential warp factors growing from the negative 
tension brane to end of the world positive tension branes. Gain in the scale due to the 
warp factor between the positive tension brane at y = and the negative tension brane 
at y = Li is lost partially by the time it reaches our universe at y = L2. It is easy to see 
that if this construction of three 3-brane system stabilizes for Li ~ L2/2 then the loss in 
the scale hierarchy is near total. It is, therefore, crucial that L1/L2 ^ 1. 

in 



In the light of this, we take up the issue of stability of the Kogan et al. model [p3 



this paper. We use the method of Goldberger and Wise[|T8[ of coupling a bulk scalar field 
to the brane system. We find that the system does stabilize, and remarkably enough, for 
L1/L2 ~ 1. We also find that this result is not very sensitive to small changes in the mass 
m of the bulk scalar field as well as small fiuctuations of the parameter k appearing in the 
warp factor as long as m ^ /c. 



2. Stabilization of the three 3-brane Model 



In this section, we will first briefiy review the Kogan et al. model and will then 
proceed with the question of modulus stabilization. As mentioned in the introduction, the 
model of Kogan et al. contains three parallel 3-branes located at Lq = 0, Li and L2. The 



^ For other as well as related proposals of stabilization of the radion modulus see |^,19-22]. 
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fifth dimension y has orbifold geometry /Z2 and Lq and L2 are orbifold fixed points. 
The action for this configuration is 

S = [d^x [ ' dyVG{2M^R - -J^f d^xV.^J -G^^ , (2.1) 

" (i) 

where Gfj,l is the induced metric on the branes and Vi are their tensions. The five dimen- 
sional metric ansatz is given by 

rfs2 _ e-^'^'^y^'n^^dx^dx'' - rldy^. (2.2) 

This metric preserves four dimensional Poincare invariance. The warp factor cr{y) is just 
a constant conformal scale factor for the induced four dimensional metric. Although it is 
constant on a given 3-brane, its numerical value is different on different 3-branes. 

For the three 3-brane model, the warp factor grows on either side of the negative 
tension brane up to the positive tension branes. The warp factor o"(y) which solves the 
equations of motion obtained from eq.(2.1) is 

a{y) = kr,{L^-\\y\-L^\) (2.3) 

when the brane tensions are Vq = — A//c, Vi = A/k and V2 = —A/k. Due to the warp 
factor, the physical mass M on the visible brane situated at y = L2 is related to the naive 
mass parameter Mq of the four- dimensional Minkowski theory by M = WMq, where, 

W = exp{-krc{2Li - L2)) (2.4) 

is the warp factor between the visible brane and the other end of the world brane. As 
mentioned earlier, to solve the hierarchy problem between the Planck scale and the SM 
scale, W should be of o(10~^^). Notice that in the three 3-brane model, if the moving 
brane located at y = Li stabilizes close to Li = L2/2 then ~ 1. In such a situation, 
we will not generate exponential scale factor in this model despite having the visible brane 
stabilized reasonably far from the other end of the world brane. Thus the ratio L1/L2 is 
as important, if not more, as the absolute stability of the three 3-brane system. 

In order to study the stability of the three 3-brane model, consider coupling a bulk 
scalar field to the system of three 3-branes. Our technique is a straightforward gener- 
alization of the GW mechanism[^. Consider the action for a bulk scalar field of mass 
m. ^ 

^Buik = ^Jd^xJ dyVGiG^^'dAXdBX - (2.5) 



where Gab is the five dimensional metric given in eq.(2.2) with a(?/) given in eq.(2.3). 



Following ref.|18], we impose the condition that the bulk scalar field satisfies certain 



boundary conditions at the location of each of the branes. The boundary potentials are 

Sl, = j d^xXiix^ - vff at y = Li, (2.6) 

Sl2 = j d'^x\2{x^ - vlf at y = L2. 

We consider only those configuration of the bulk scalar which solve the equations 
of motion subject to the condition that the boundary potentials (at y = 0, Li and L2) 
are minimised. This essentially amounts to neglecting dynamics of x along the directions 
tangential to any of the 3-branes. This assumption is reasonable because we are looking 
only at the stability of the three 3-brane system at the moment and are not studying the 
phenomenologcal consequence of possible coupling of the bulk scalar field x to matter fields 
living on the branes. It, therefore, suffices to concentrate on the equation of motion of x 
only in y direction. Equation of motion for x along y is 

dlx - Acy'{y)dyX - m'rlx = (2.7) 

where cr'{y) = da{y)/dy. It is straightforward to find the solution to this equation and it 
is given by 

X{y) = exp {2a'{y)y)[Aex.p {a'{y)uy) + Sexp {-a'{y)uy)] (2.8) 

where, v = ^4 + (mrj Va'2(y), (i, B) = (A, S) for < y < Li and (i, B) = (C, D) for 
Li < y < L2. It is worth mentioning here that the formal solution (2.8) is a function of 
cr'(j/) but within any given range of values of y, e.g. < y < Li, a'{y) is independent of 
y. The coefficients A, B, C and D are determined by demanding that x minimizes the 
boundary potential. This gives, 

(2 9) 

where, for notational convenience, we have made a change of variables from Li and L2 to 
X = exp {—krcLi) and Y = exp {—krcL2). Substituting this solution into the action and 
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integrating out y gives a four dimensional effective potential for X and Y. Writing this 
effective potential in terms of X and a new variable R = Y/X, we have, 



An important thing to notice at this point is that, for arbitrary positive values of i^, 
the potential (2.10) grows as X — > 1 or as i? — > 1 as long as vi ^ vq and vi ^ V2. These 
two limits correspond to the negative tension brane approaching the positive tension brane 
at J/ = and sX y = L2 respectively. In fact, the potential has the following singularity in 



these limits: V{X,R) ~ (1 - X^'^)-! > as X ^ 1 and V{X,R) ~ (1 - i?^'^)-! > as 



i? — > 1. This implies that the negative tension brane experiences repulsive forces exerted 
on it by both the positive tension end of the world branes and thus the three 3-brane 
model cannot reduce to the two 3-brane model. In other words, for a fixed value of L2, L\ 
lies strictly inside the interval (0,L2). However, there is an interesting caveat in this and 
that has to do with the choice of vev's of the bulk scalar field on the 3-branes. If vi = vq 
and/or vi = V2, then the leading singularity in the potential is removed and the subleading 
terms in the potential are attractive. In other words, for vi = ^2, the potential (2.10) is 
such that for R close to 1, V is attractive and assumes a finite value for i? = 1. In this 
case, therefore, the two 3-brane limit is a stable one. The situation is similar, though not 
identical, for vi — vq and X — > 1. Prom here onwards we will work with the case when vq, 
vi and V2 take different numerical values. We still need to find out whether it is possible 
to stabilize L2. An equally important question relates to the magnitude of L2 — Li. This 
number is important to get the correct metric scale factor on the visible brane. 

We will first show how the three 3-brane system is stabilized by the bulk scalar field. 
In the process we will also be able to determine L2 — Li. Change of variables from X, 
y to X and R has simplified the form of the potential (2.10) to a considerable degree, 
yet, minimization remains a complicated task. The first point to notice here is that the 
potential is bilinear in Vi and hence the results depend only on the ratios 

ro = — , r2 = — . (2.11) 

Vi Vi 

Furthermore, this choice of variables is such that the Hessian matrix is automatically 
diagonal on the extremal locus. Extremizing V with respect to R, we get 




(2.10) 




(2.12) 
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Fig. 1: As the ratio r2 of the vev of x the visible brane to that on the 
negative tension brane approaches 1, location of these two branes coincides. 
However, for < 1 we have two possible values for r2 as can be seen in this 
graph. 



An analytic inversion of this equation is not possible on account of its transcedental 
nature. Numerical solutions are straightforward though and we present the results in Fig.l. 
For this purpose, instead of working with z/, it is instructive to consider 

2 

T77 

e^.-2^— , (2.13) 



with the approximate equality holding for a light bulk scalar field. 

Although the solutions of eq.(2.12) exist for a wide range of e, we are interested in 
small e's and shall limit ourselves to this regime. In Fig.l, the left and right branches 
correspond to and respectively. At the cusp (r2 = 1), the solutions move into 
unphysical domains. Looking at the figure, we surmise easily that Rk, 1 demands r2 ~ 1. 
This is particularly pronounced for small e. 

While the solutions eq.(2.12) represent minima in the i?-direction, it remains to be 
seen whether simultaneous minima in the X-direction exist. We answer this question next. 
Extremizing with respect to X, we obtain 



(2.14) 
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where 



B=(l-i?2^) (zy2 + 4zy-12)-8z/ + 16zyi?''-V2 + 16ri-4r|(z/ + 2) + . (2.15) 

In each of eq.(2.14), r2 assumes both of the values r^. We thus have four possible branches 
in the solution space. Numerically though, for given R and X, the four solutions for 
ro roughly split into two pairs (the two r^s on the one hand and the two TqS on the 
other), with the intra-pair splitting considerably smaller than the inter-pair one. More 
interestingly, for a given e, 

• not the entire curve of Fig. 1 is admissible when confronted with eq.(2.14). Rather, only 
r2 ~ 1 (and hence ~ 1) can lead to a stable solution. For r2 significantly different 
from unity, the quantity B in eq.(2.15)becomes negative resulting in unphysical values 
for X; 

• for ro = Tfj", d'^V/dX^ < irrespective of the choice for r2- Thus, each of these 
branches correspond to a sequence of saddle points. 
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Fig. 2: (a) For small values of e, the warp factor is not very sensitive to 
the vev ratio r2. The dependence on tq is rather pronounced though. Warp 
factor of 10~^^ can be obtained with both the ratios ro and r2 of o(l). (6) 
For small e, a warp factor of 10~^^ can be obtained with both the ratios ro 
and r2 of o(l). However, for larger values of e, the ratio ro is required to be 
quite large. The numbers in the parentheses refer to (e, N) where is the 
factor by which the y-axis has been rescaled for each curve. 
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In Fig. 2(a), we exhibit the relation between vq and r2 that must be satisfied to obtain a 
particular warp- factor for a given e. As is evident, is a very sensitive function of vq. 
This is not very unexpected as, for small e, eq.(2.14) essentially gives ro to be polynomial 
function of \nX (and hence Inl^). This exponential dependence, a feature we share with 
the GW solution to the original RS model, could of course be termed a weakness of such 
stabilization schemes. 

The dependence on e — see Fig. 2(6) — is even more pronounced. This, again, is not 
unexpected as e ~ o(l) implies m ~ A;. For such large masses of a field propagating in the 
anti-de Sitter bulk, its wavefunction decays very fast. Consequently, its classical values at 
the two end of the world branes would be widely different. 



3. Discussion 

Brane world universe is one of the promising proposals for exploring the physics beyond 
the Standard model. There are several variants of the original proposal of Randall and 
Sundrum0. Many of them are invoked to aviod possible contradications with our existing 
knowledge of the Standard model physics and the standard big bang cosmology. The 
proposal of Kogan et al. is along the same lines. We have shown in this paper that the 



three 3-brane model can be stabilized by coupling the system to a bulk scalar field and 
using the Goldberger-Wise formalism [|18| . 



We find that stability of the three 3-brane model is a more delicate problem than that 
for the Randall-Sundrum model. This is because gain in the scale due to the warp factor in 
the Randall-Sundrum model is partly offset by the moving brane. It is therefore not only 
necessary to have overall stabilization of the three 3-brane system but it is also necessary 
to have the ratio L1/L2 close to 1. We show that the generalization of the Goldberger-Wise 
mechanism to this model stabilizes the radion modulus in such a way that, for relatively 
small mass of the bulk scalar field or, equivalently, small e, L1/L2 ~ 1. Consequently, it 
is possible to generate a scale hierarchy between the Planck scale and the TeV scale even 
without finetuning the parameters. 

As we have discussed earlier, the modulus ceases to be stabilized in the event that the 
vevs on two (adjacent) branes are exactly equal. Such a degeneracy might seem natural, 
at least in the classical limit. However, it should be realized that quantum corrections 
are likely to disturb such an equality. This is particularly so since the vev on a brane is 
affected, at one loop level and higher, by the particle spectrum on the brane and their 
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coupling to the bulk scalar field. Unless the spectrum is exactly alike on two adjacent 
branes, and unless their coupling to x are exactly the same, it seems unlikely that the two 
vevs could be equal. And since it is only in this exactly equal vevs limit, that the potential 
becomes attractive, the point is, perhaps, a moot one. 

Acknowledgments: We would like to thank A. Sen for discussions. 
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